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\ We present an analysis of X-ray high quality grating spectra of the Seyfert 1 

galaxy NGC 5548 using; archival Chandra-'H.'ETGS and LETGS observations for 
a total exposure time of 800 ks. The continuum emission (between 0.2-8 keV) 
q \ is well represented by a power law (T = 1.6) plus a black-body component (kT 

| =0.1 keV). We find that the well known X-ray warm absorber in this source 

ch \ consists of two different outflow velocity systems. One absorbing system has a 

velocity of -1110 ± 150 km s _1 and the other of -490 ± 150 km s _1 . Recognizing 
the presence of these kinematically distinct components allows each system to be 
| fitted independently, each with two absorption components with different ioniza- 

tion levels. The high velocity system consists of a component with temperature 
of 2.7 ± 0.6 x 10 6 K, logU = 1.23 and another component with temperature of 
5.8 ± 1.0 x 10 5 K, logU = 0.67. The high- velocity, high-ionization component 
^ ; produces absorption by charge states Fexxi-xxiv, while the high- velocity, low- 

ionization component produces absorption by Neix-x, Fexvil-xx, Ovil-vill. 
The low-velocity system required also two absorbing components, one with tem- 
O perature of 5.8 ± 0.8 x 10 K, logU = 0.67, producing absorption by Neix-x, 

Fexvil-xx, O VII- vill. The other with lower temperature (3.5 ± 0.35 x 10 4 K), 
and lower ionization (log £7 = —0.49); producing absorption by Ovi-vn and the 
Fe vil-xil M-shell UTA. Once these components are considered, the data do not 
5_i ■ require any further absorbers. In particular, a model consisting of a continuous 

radial range of ionization structures (as suggested by a previous analysis) is not 
required. 



1 Instituto de Astronomia, Universidad National Autonoma de Mexico, Apartado Postal 70-264, 04510 
Mexico DF, Mexico. 

2 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street. Cambridge MA 02138, USA. 

3 Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA. 

4 OAR-INAF, Via Frascati, 33, 00040, Monteporzio Catonc, RM (Italy) 

5 IESL, Foundation for Research and Technology, 711 10, Hcraklion, Crete (Greece) 



-2 - 



The two absorbing components in each velocity system are in pressure equi- 
librium with each other. This suggests that each velocity system consists of a 
multi-phase medium. This is the first time that different outflow velocity systems 
have been modelled independently in the X-ray band for this source. The kine- 
matic components and column densities found from the X-rays are in agreement 
with the main kinematic components found in the UV absorber. This supports 
the idea that the UV and X-ray absorbing gas is part of the same phenomenon. 
NGC 5548 can now be seen to fit in a pattern established for other warm ab- 
sorbers: 2 or 3 discrete phases in pressure equilibrium. There are no remaining 
cases of a well studied warm absorber in which a model consisting of a multi-phase 
medium is not viable. 

Subject headings: galaxies: absorption lines - galaxies: Seyferts - galaxies: active 
- galaxies: X-ray 

1. Introduction 

In ~ 50% Seyferts 1 and quasars there is an interesting phenomenon, the ionized ab- 
sorber or warm absorber (WA) (Reynolds 1997; Crenshaw et al. 2003a; Piconcelli et al. 
2005). This phenomenon was observed in the X-ray band for the first time in the spectrum 
of the quasar MR2251-178 by Halpern (1984) with the Einstein Observatory. More recently, 
high resolution observations by Chandra and XMM-Newton have revealed the presence of 
blueshifted absorption lines (~v out ~ —500 to —2300 km s -1 with respect to the AGN host rest 
frame; Collinge et al. 2001; Kaspi et al. 2002), produced by warm ionized gas (T~ 10 4 — 10 6 
K; Krongold et al. 2003). These reveal the nature of the WA as an outflowing wind. The 
absorbing gas usually shows several components with different states of ionization, as it is 
demonstrated, for instance, by the different transitions of Oxygen (O VII, O vm), or the 
different transitions of Iron (FeVll-xil and Fexvil-xxil). Transitions by Civ, Nv, and 
Ovi are observed in both the X-ray and UV spectra of these sources (with similar outflow 
velocities, suggesting a connection between the narrow absorption line systems in the UV 
and the WA (Crenshaw et al. 1999, Kaastra et al. 2000; Kaspi et al. 2002; Krongold et al. 
2003; Steenbrugge et al. 2005). 

The possibility that the absorbing components with different ionization states are in 
pressure equilibrium with each other, has been suggested by the analysis of several Seyfert 
1 galaxies: NCG 3783 (Krongold et al. 2003, Netzer et al. 2003), NCG 985 (Krongold et al. 
2005; 2009), NGC 4051 (Krongold et al. 2007), Mrk 279 (Fields et al. 2007). These results 
strongly suggests that the absorbers consist of a discrete multi-phase medium. Is important 



-3 - 



point out that these studies were carried out by fitting observations using photoionization 
models. A detailed model for the WA consisting of a thermal multiphase wind, with a density- 
size spectrum has been presented by Chelouche & Netzer (2005); Chelouche (2008). Another 
representation of the absorber (which includes calculations of line radiative transfer, but does 
not fit directly the data) may be a constant pressure medium with a radially stratified density 
distribution (Rozahska et al. 2006). 

1.1. The Warm Absorber in NGC 5548 

NGC 5548 (classified as a Seyfert 1 with Broad Emission Lines, Full Width Half Maxi- 
mum (FWHM^) = 5610 km s _1 ; Turner et al. 1999) is at a redshift z = 0.017175±0.000023 
(De Vaucouleurs, 1991; the determination of z was performed with the detection of an HI 
emission line with negligible error in the velocity, Heckman et al. 1978), and has a luminosity 
of L 2 _ 6KeV = (2.0 - 2.8) x 10 43 erg s^ 1 (Branduardi-Raymont, 1986; Nicastro et al. 2000; 
present work). 

This galaxy has been studied many times in the optical, UV, and X-ray bands. The 
optical variability of the source allows detailed reverberation mapping studies of the Broad 
Emission Line Region (BELR; e.g. Peterson, 1993), yielding a central black hole mass of 
6.7 ± 2.6 x 10 7 M (Peterson 2004). This is ~ 33 times larger than, for instance, the 
NGC 4051 black hole mass (2xl0 6 M , Peterson et al. 2004)0, for which the WA has been 
extensively studied and has been described as an accretion disk wind (Krongold et al. 2007). 
The BELR radius calculated using the H/3 emission line for NGC 5548 is 21.2±2.4 1-d (Kaspi 
et al. 2000) and the He II BELR is 16.2 1-d (Bottorff 2002). This is ~ 3 and - 5 times 
larger, respectively, than the NGC 4051 BELR radius. 

The X-ray flux is variable (Branduardi-Raymont, 1989; Nandra et al. 1991; Chiang et 
al. 2000; Markowitz et al. 2003). The long-term X-ray and optical continuum light curves 
are highly variable and are correlated (Uttley et al. 2003). It has been suggested that these 
variations are due to thermal instabilities in the inner disk (Treves et al. 1988). The short- 
term X-ray variability may be more concentrated toward the center of the accretion flow 
(Uttley et al). 

Low resolution X-ray studies (e.g. Nicastro et al. 2000, using BeppoSAX satellite ob- 



1 However, note that Marconi et al. (2008) suggested that the black hole mass of Narrow Line Seyfert 1 
galaxies (like NGC 4051) is underestimated by a factor ^5, when radiation pressure is taken into account in 
the virial theorem. 
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servations) showed the presence of an ionized absorber, as the data presented transitions 
consistent with bound- free absorption by H-like and He-like ions of O. The first high resolu- 
tion spectrum of NGC 5548, obtained with the Chandra Low Energy Transmission Grating 
Spectrometer (LETGS), confirmed the presence of the ionized absorber with many absorp- 
tion lines from transitions such as C VI, Nvn, O VII- vill, and Neix-x (Kaastra et al. 2000). 
In the UV band, spectra obtained with the Goddard High-Resolution Spectrograph (GHRS) 
on the Hubble Space Telescope (HST) presents multiple velocity components, with intrinsic 
absorption lines such as Ly a , CivA1549, and NvA1242 (Mathur et al. 1999, Crenshaw et 
al. 2003b). Five different kinematic blueshift components are observed with velocities cor- 
responding to -166, -336, -530, -667, and -1041 km s _1 (components 5 to 1 in Crenshaw et 
al. 2003b). Far Ultraviolet Spectroscopic (FUSE) spectra also show intrinsic Ovi and Lyg 
absorption (Brotherton et al. 2002). 

NGC 5548 is the first Seyfert galaxy for which a possible connection between the X- 
ray and UV absorbers was suggested (Mathur et al. 1995) given the ionization state of 
the medium (though a connection has already been suggested for the quasar 3C 351 by 
Mathur et al. 1994). However, the exact nature of this connection in NGC 5548 has been 
controversial. Brotherton et al. (2002) disputed this conclusion by the discrepancies in 
the Ovi ionic column densities inferred in the two bands. Arav et al. (2003) then found 
consistency between these columns, once the effects of velocity-dependent covering fraction in 
the UV band are considered. However using simultaneous Chandra-LFiT GS and HST-STIS0 
observations, Crenshaw et al. (2003)b report that the ratio j^7^3 , in the UV component 
1, is 2.5, which suggests a gas with low ionization state, implying low O vil and O vm 
column densities. Crenshaw et al., then suggested no relation between the UV and X-ray 
absorbers. Using Chandra High Energy Transmission Grating Spectrometer (HETGS) and 
LETGS data, Steenbrugge et al. (2005) proposed again a link between that X-ray and UV 
absorbers, given the similar kinematics seen in the absorption lines in the two bands. A 
new study based on 2002 UV spectra suggested that at least a part of UV component 3 
should produce absorption in X-rays range (Crenshaw et al. 2009). While it is clear that 
some connection exists, the exact relation between the UV and X-ray absorbers remains 
uncertain. 

NGC 5548 is perhaps an exception to the idea that the WA is a multi-phase medium 
possibly in pressure balance (JQ). A continuous radial range of ionization states was suggested 
for the structure of the WA in this source (Steenbrugge et al. 2005). It should be noted that 
NGC 5548 has an unusually low Galactic column density (1.6 xlO 20 cm -2 , Murphy et al. 



2 (Hubble Space Telescope- Space Telescope Imaging Spectrograph) 



-5 - 



1996), so the X-ray spectrum remains less absorbed at soft X-ray wavelengths than in most 
AGNs. This broader wavelength range exposes extra atomic transitions that can provide 
additional constraints on WA conditions (for instance NGC 3783, has a very well studied 
absorber, but has a substantially higher Galactic column density than NGC 5548, 9.5 xlO 20 
cm -2 , Murphy et al. 1996). Thus, it is worthwhile to test the multi-phase scenario for the 
absorber in this object. 

For these reasons, an analysis of all the seven Chandra datasets seems in order. We 
re-analyzed archival Chandra LETGS and HETGS high resolution spectra of NGC 5548. 
Despite previous results that suggest that only a continuous radial flow could model the 
absorber (Steenbrugge et al. 2005), we show that its physical structure is also consistent 
with a multi-phase medium. We contrast our results with previous studies of this source. In 
section 2 we describe the observations and data reduction. In section 3 we describe our data 
modeling. In § 4 we discuss our results, and in § 5 their implications. 



2. Data Processing 

NGC 5548 was observed seven times with the Chandra X-ray Observatory (Weisskopf et 
al. 2000) between 1999 and 2005 (see Tabled]) . Two observations were carried out with the 
HETGS (Canizares et al. 2000) with the light dispersed on to the Advanced CCD Imaging 
Spectrometer (ACIS; Garmire et al. 2003). The remaining five observations were obtained 
with the LETGS (Brinkman et al. 2000) and the High Resolution Camara (HRC; Murray et 
al. 2000). We retrieved the primary and secondary data products of the observations from 
the Chandra data archival The data reduction was performed with the Chandra Interactive 
Analysis of Observations software (CIAO v.3.4 0; Fruscione et al. 2002). We followed the 
data analysis threads provided by the Chandra X-ray CenteJE Negative and positive first 
order spectra were extracted and response matrices created for each observation. In this 
paper, we explore the time-averaged properties of the spectra, focusing on the best possible 
determination of the different ionization and velocity components of the well-known ionized 
absorber in this AGN. Thus, we co-added the spectra obtained with each instrument to 
get the maximum signal to noise ratio (S/N). This gave a total net exposure time of 236 
ks for the MEG and HEG spectra, and 564 ks for the LEG spectrum. The source shows 
flux differences among observations. The change in flux between the 2 HETGS observations 
is pal. 7. The change in flux among the LETGS data is a factor of 6 during the six years 



3 http://cxc. harvard.edu/cda/ 

4 http://cxc. harvard. edu/ciao3.4/indcx. html 
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of observations (see FigJT]). The change in flux in short time scales between consecutive 
(LETGS) observations is wl.2-1.6 (see Tabled]), (for observations separated by 3 days). A 
time-evolving analysis is postponed to a forthcoming paper (Andrade- Velazquez et al. in 
preparation). However, we discuss the possible effects that these variations in our analysis 
in Section 13.2.21 

We have analyzed the data in the range between (1.57 to 10.81) A ([1.15 to 7.90] keV) 
for the HEG, (2.46 to 24.58) A ([0.50 to 5.04] keV) for MEG, and (8.85 to 49.15) A ([0.25 
to 1.40] keV) for LEG. All the data were binned by a factor of 2, corresponding to 0.005 A, 
0.01 A, and 0.025 A per bin for the HEG, MEG and LEG, respectively. 

3. Modeling 

We have used the CIAO software Sherp80 (Freeman et al. 2001) to perform spectral fits 
to the data. In all the models described below, the attenuation due to the Galactic absorption 
was accounted for with a neutral absorber (see § 1). We first focus on the HETGS (HEG 
and MEG) data. Then, we compare the final HETGS model with the LETGS data, and find 
a good fit at both short and long wavelengths (see below). 

3.1. Fitting the HETG data. 

3.1.1. MODEL A: No Absorption Components. 

Initially, we fit the data with a power law, following previous analysis (Reynolds et 
al. 1997: Advanced Satellite for Cosmology and Astrophysics -ASCA; Nicastro et al. 2000 
BeppoSAX) and a black-body component (Kaastra & Barr 1989, Steenbrugge et al. 2005), see 
Table [21 Model A. We further added four Gaussians to model the evident emission lines (see 
Table 0Q) . 

Three emission lines are resolved, with FWHM between 500-700 km s _1 (MEG has an 
FWHM resolution of ±2 155 km s" 1 at 19 A). These lines are in the soft band (consistent 
with NelxA13.69, O VllA 22.10, and NviA 23.27). The short wavelength « 5 A) region of 

Tittp://cxc. harvard.edu/ciao/threads/indcx. html 
6 http:/ /cxc. harvard.edu/sherpa 

7 This Table lists eight emission lines, as a detailed analysis of the spectra §3.4 revealed the presence of 
another four emission lines 
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the spectra shows the presence of the fourth emission line. This corresponds to the commonly 
found low-ionization (< Fexvil) Fe K a line (6.4 keV, Nandra et al. 2007), and is located 
at 1.937 A (without any shift from the rest frame of the host galaxy). Yaqoob et al. (2001) 
report an EW of 133 eV and a width of 4515^645 kms' 1 and suggest a possible origin in 
the outer Broad Line Region (BLR). We report a full width half maximum (FWHM) of 
2644±2469 km s -1 consistent with the value reported by Yaqoob et al. We find a slightly 
lower EW than this authors, EW~ 92.2±43.9 eV, however, a detailed analysis of this line is 
beyond the scope of this paper. 

Though model A was statistically acceptable (x 2 /d.o.f.= 3813/4112), strong negative 
residuals between (5-24.5) A are present, consistent with absorption features by ionized 
atoms (Fe xvii-xxii, Ovii-vin, Neix-x, Sixm-xiv, Nvi-vii, and Mgxi-xn) due to the 
well-known WA in this source (see Fig [5]). 

3.2. Ionized Absorber : fitting with PHASE 

To model the features produced by the ionized absorber, we used the code PHASE 
(Krongold et al. 2003). The code has 4 free parameters per absorbing component, namely 
1) the ionization parameter at the illuminated face of the absorber U = ■> where Q(H) 

is the rate of H ionizing photons, r is the distance to the source, uh is the hydrogen number 
density and c is the speed of light; 2) the hydrogen equivalent column density Nh, 3) the 
outflow velocity of the absorbing material V out and 4) the micro-turbulent velocity Y tU rb 
of the material. Usually, the micro-turbulent velocity is not left free to vary, because the 
different transitions due to ionized gas are heavily blended, or because the different velocity 
components are also blended and cannot be resolved (e.g. Krongold et al. 2003, 2005 for 
NGC 3783). In the case of NGC 5548, with the HETGS data resolution, it is possible to 
distinguish two velocity components (see below). So, despite the fact that the individual 
absorption lines cannot be resolved, we have left the micro-turbulent velocity free to vary 
(at least for the strongest absorbing components, see below). PHASE has the advantage 
of producing a self-consistent model for each absorbing component, because the code starts 
without any prior constraint on the column density or population fraction of any ion. We 
have assumed solar elemental abundances (Grevesse & Noels 1993). So, given the intrinsic 
Spectral Energy Distribution (SED, the approximation of this input parameter is explained 
in the next subsection) of the source, the column density, and the ionization parameter of 
the absorbing media, the ionization balance is fixed. 

In all the models described below, we have assumed that the continuum source (including 
the soft excess) is fully covered by the ionized absorbing gas while the emission lines are 
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produced outside the region where the absorbing material is located. 

3.2.1. Spectral Energy Distribution 

The X-ray to radio SED of NGC 5548, used to calculate the ionization balance of the 
gas, was defined in the following way: Between the radio and the near UV ranges, we used 
(non simultaneous) data from NED (NASA/IPAC Extragalactic Databas^l, references are 
attached in Appendix In the range between 100 fxm and 912 A, we performing a linear 
fit to the over 100 references of flux measurements found in NED. For the region between 
(0.25-7.90) keV ( [1.57 - 49.15] A), we used the power law (T) as well as the black-body 
contribution, as inferred from the fits to the continuum presented in Table HI For energies 
> 8 keV, we extrapolated the derived T in our fits up to 130 keV, where a cut-off is observed 
(Nicastro et al. 2000). Given that the UV SED shape has only a second order effect on 
the ionization balance of the X-ray species, as these charge states have ionization potentials 
larger than 0.1 keV (see Netzer 1996, Steenbrugge et al. 2003), we have connected the 
unobservable region of the SED between 912 A and 0.25 keV with a simple power law (Haro- 
Corzo et al. 2007). The radio region (at 21 cm) was connected with the infrared (100 /im) 
also with a simple power law. Figure [3] shows the final SED used in this work. 

3.2.2. Modelling Averaged Spectra 

An important caveat of averaging the different datasets is that the the physical prop- 
erties of the absorbers might have changed among observations, either because of intrinsec 
changes or because of changes induced by the variations in the flux level. We note, that 
Kaastra et al. 2002 and Steenbrugge et al. 2005 report similar values for the Hydrogen 
equivalent column density (within the uncertainties) for the 1999 and 2002 observations (see 
Table [T|). Then, we do not expect a significative change in the equivalent Hydrogen col- 
umn density among the seven observations. The same is true for the outflow velocity of the 
absorbing gas. 

However, the flux level changes between the observations used in this work. Assuming 
the gas is in photoionization equilibrium, this would drive changes by the same factors 
observed in the ionization parameter (U) of the gas (producing changes in the gas opacity; 
Nicastro et al. 1999). There is a flux variation by a factor of 1.7 between the HETG 



http : / /nedwww .ipac.caltech.edu/ 
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observations, thus, the expected variation between these two observations will be Delta 
(logU) = 0.23. For the LETG observations, there is a large change in flux (by a factor of 6 
during 6 years; see Fig. [T]) between the first and last two observations. In photoionization 
equilibrium this would produce a change Delta(logU) =0.77 in the ionization parameter. 
Thus, our determination of this parameter in the co-added data will represent a weighted 
average over the different observations. We note that for the LETG, the 1999 and 2002 
observations have the higher S/N and similar fluxes (different by a factor ~ 2.19, and also 
similar to the HETG observations). Then, the best fit value of the absorbing gas ionization 
parameter on the LETG co-added spectrum will be highly dominated by these observations 
(the inclusion of the 2005 fainter observations contribute to increase the S/N by ~ 10%). 



3.3. Ionized Absorber 

3.3.1. Model B: One Ionized Absorber Component. 

We made a new model, B, by including a single ionized absorbing component to model 
A. Table [5] shows the best fit parameters values for this model. The outflow velocity is 
~ -740±150 km s _1 and could correspond to the UV component 2, which has a value 
~ -667 km s _1 . Given the ionization parameter (logU^) and column density (logN^s) of 
this component (see Table [2]), this model fits the absorption lines produced by ions such as 
Fe xvii- xxii, Neix-x, Sixill, Mgxi-xil, and O vm (see Fig. H]). The fit gives a significant 
improvement over model A (x 2 = 3053 for 4108 d.o.f.; A% 2 = 760). An F-test gives a higher 
than 99.99% confidence level for the presence of the absorber (see Table [2]). The absorbing 
component can reproduce ~ 20 absorption features (including blends) in the wavelength 
region (5 - 20) A. However, significant residuals remain at 6.18 A, which are probably due 
to a Si XIV transition (A6.182), and at 10.62 A and 11.01 A, probably corresponding to 
transitions by Fexxiv(A10.619) and Fexxill(A11.019), see Figure [5) indicating the possible 
presence of another absorbing component of gas with higher ionization. Strong residuals 
found close to 16.5 A, could correspond to inner transitions of Feix, further suggesting the 
possible presence of another absorber with lower ionization state. 

In addition, a detailed inspection of the spectrum near the absorption lines fit by this 
model (Mgxi, Neix, Fexvil, and Sixiv) reveals a complex structure in the line profile, 
with additional absorption with an outflow velocity of ~ -1040 km s _1 (see Fig. [6]). These 
results, implying additional components and complex velocity structure, are consistent with 
the analysis of Steenbrugge et al. (2005). We will model these additional components in the 
following sections. 
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3.3.2. Model C: Two Ionized Absorber Components. 

Model C includes two absorbing components, whose parameters are again free to vary. 
This model gives a x 2 — 2987 for 4104 d.o.f., significantly improving over model B (A% 2 = 66, 
an F-test, gives a confidence of larger than 99.99% for the addition of these parameters). 
The ionization parameters (logUci and logU<7 2 ) are different by almost a factor of 3 (see 
Table E]). 

The component with lower ionization (with logUci slightly lower than logUci) models 
the transitions already fitted in model B: Si xm-xiv, Ne ix-x, Mg xi-xil, Fe xvil-xxil, and O 
vm (see Fig. [7j). The higher ionization component (with logUci ~ 1. 4 x logUci) contributes 
weakly to fit these ions, but fits the transitions due to Fexxm-xxiv (not fitted in model B). 
The two components have different outflow velocities ( V ou tci — -1120 km s _1 and V ou tc2 — 
-450 km s _1 ) . We note that while model B partially fits the Si xiv line, the higher ionization 
component in model C (logUci), which has also higher outflow velocity, improves the fit to 
the Si xiv line, i. e. this component fits the residuals that model B left around 6.18 A (see 
Fig. [7]). While model C does reproduce, then, the lines produced by Si xiv , Fexxi-xxiv, it 
cannot model the complex absorption profile observed also for other lines. 

3.3.3. Model D: Three Ionized Absorber Components. 

Motivated by the line profile residuals, we built a new model, D, by adding a third 
absorbing component, whose parameters are again free to vary, as are those of the other 
two components. Model D gives a \ 2 = 2946 for 4100 d.o.f., improving significantly model 
C (Ax 2 = 41, F-test confidence level of 99.99%). All the three absorbing components have 
different ionization level. Two of the three absorbing components have similar values of 
logU to those found in model C. The third component has a much lower ionization level 
(by at least a factor ~ 14, see Table E]), and fits the Nvi-vil, O VII, Mgix-xi lines, as 
well as the FeVll-xil M-shell Unresolved Transitions Array [UTA] (see Fig. [8]). The three 
absorbing components group in two different outflow velocity systems. One system, with 
outflow velocity in the range between -450±150 and -590 ±150 km s _1 , contains the systems 
with intermediate ( logU^) and low (logU^) ionization level. The other velocity system 
(V OU ( = -1116 ± 150 km s _1 ) is formed by the high ionization component. 

We note that model D still leaves some residuals around 16.5 A, where an Fe IX transition 
lies. The reason for this discrepancy may be due to our use of the abbreviated data to model 
the UTA (see Krongold et al. 2003 and Behar et. al 2001 for further details). 
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3.3.4- Model E: Four Ionized Absorber Components. 

Model D does not solve the problem of the line profile for several absorption lines with 
intermediate ionization state, in particular the Fexvn and Nex transitions. Therefore, 
model E was built by adding a fourth absorbing component to model D. This case gives a 
X 2 = 2904 for 4096 d.o.f. (Ax 2 = 42), which is statistically better than that of model D (an 
F-test shows that the fourth absorber is required at the 99.99% confidence level). The best- 
fit parameters value for each absorbing component are listed in Table [6j and the global fit to 
the HETGS data is shown in Figures M (for MEG data) and M (for HEG data) . Model E 
again consists of two outflow velocity systems (hereafter High Velocity or HV for the system 
with V out = -1110 ± 150 km s" 1 , and Low Velocity or LV with V out = -490 ± 150 km s" 1 
), each having two different ionization components^. The ionization parameter of the high 
ionization phase in the LV system (hereafter LV-HIP) is logU=0.67 with a column density 
of \ogN H = 21.25. The low ionization phase (LV-LIP) has \ogN H = 20.74 and logU=-0.49. 
The features fitted by this system are the same ones fitted in model D by the components 
similar to that model. The contribution of this velocity component to the global fit is shown 
in Figure [TTJ 

The HV system has two absorbing components. One of the absorbing components of 
this system has a high ionization level (logU=0.671, hereafter called HV-HIP), and a column 
density of logN#=21.03, and the other absorber has an even higher ionization level (here- 
after called HV super-high ionization phase or HV-SHIP) with logU=1.23 and logN#=21.73. 
Figure [12] shows the contribution of the HV system to the global fit of the spectrum. The 
HV-SHIP fits the same features model D does. 

The HV-HIP of this system fits the same absorption lines as the LV-HIP (given that the 
ionization parameters of these two components are quite similar). The difference between 
these two components lies mainly in their outflow velocities, finally improving the fit to the 
line profile (as seen in Fig. [13]) . We note that this result is obtained without any constraint 
in the main parameters of the absorbers, which further confirms the presence of the two 
velocity systems. Separating the two components was possible given the excellent resolution 
of the HETGS, ^155 km s" 1 (19 A) for MEG0 



9 The velocities for the HV and LV systems were calculated taking the average of the two components 
forming each system (see Table |6|) 

10 The resolution here and anywhere else is given with respect to the factor of 2 binning used in this 
analysis. 
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3.4. Emission Features in the Spectra of NGC 5548 

In the previous models we fitted four strong emission lines (the Fe K a line, and the Ne IX, 
O VII, and N VI forbidden transitions). However, in the MEG data, there are four additional 
emission lines that can be detected with a significance level > 2a once the absorbers are 
properly modelled. These lines correspond to the resonant and intercombination transitions 
of O VII, the resonant transition of N VI and Si xill line. We produced a new model including 
these lines (model F). The best fit parameters for the eight absorption lines included to fit 
the HETG data are given in Table [3j Figure [9] already includes these features. 

Our analysis does not require the inclusion of broad emission features (as was reported 
by Steenbrugge et al. 2005). The forbidden (z), resonant (w), and intercombination (x+y) 
lines of He-like ions (Gabriel & Jordan 1969) can be used to derive the physical conditions of 
photoionized plasma (Porquet & Dubau 2000). In particular, the ratio between the resonant 
and intercombination lines serve as a diagnostic of the electronic density n e (R(n e ) = ^^). 
Using the R value of the O VII lines (R(n e )= 2.23±1.6), and following Porquet & Dubau 
(their Figure 8), we find that the gas must be photoionized and must have a density n e = 
2.0 ± 0.8 x 10 10 cm~ 3 . Kaastra et al. (2000) report an upper limit for n e < 7.0xl0 10 cm~ 3 , 
with a R(n e ) = 2.22±3.45, which is in agreement with our result. 



3.5. Extrapolation to LETGS 

The final model (F) for the ionized absorber in the MEG and HEG spectra of NGC 
5548 can be further tested by extrapolating to the LETGS data (the LETGS resolution is a 
factor ~ 2.5 worse than the MEG resolution at 19 A). In doing this, we implicitly assume no 
significant opacity variations between the co-added HETGS and LETGS spectra that may 
arise from flux variability (Nicastro et al. 1999, Krongold et al. 2005b, 2007). Given that 
the average flux level in these two spectra differs by only 30% (Steenbrugge et al. 2005), this 
is a reasonable assumption 0. 

To fit the LEG spectrum, we used model F, with the 4 absorbing components' parame- 
ters fixed to the values obtained for the HETGS data (hereafter model G). We also used the 
same continuum components (a power law and a black-body) , but we allow these parameters 
to vary freely in the fit, as the continuum level and spectral shape might be different (in 



n This is not the case between individual observations, where opacity changes might be present. See 
Andrade- Velazquez et al. (in preparation) for a full analysis of variability of the absorber to continuum 
changes. 
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addition, the LETGS is more sensitive to low energy photons and so these data allows for 
a better determination of the black-body temperature). In model G we also let free to vary 
the parameters for the emission lines detected in the HETGs spectra (except for the Fe ka 
and Si xm lines that lie outside the spectral range of the LETG). We further added three 
additional emission lines, (located in wavelengths larger than 25 A, the lower bound of the 
MEG data), corresponding to NvA29.502, NviA29.029, and CvA41.429). 

An acceptable fit is obtained to the LETGS spectrum with this simple method (x 2 / d.o.f. = 
2013/1618), implying that the HETGS spectra alone can constrain in a reliable way the prop- 
erties of the ionized absorber. The best-fit continuum parameters for the LETGS spectrum 
are collated in Table HI Table [3] presents the parameters for the emission lines. 

The best-fit model over the data is presented in Figure (see also Table [2]). This 
figure show that most absorption lines observed in the HETG data are well reproduced by 
the model in the LETGS spectrum. It also shows that model G fits well more than 12 low 
ionization absorption lines in the region (26-49) A. These lines correspond to transitions 
of ions such as Cv, Nv-vi, and Sxn. In the UTA region there are discrepancies between 
data and model, but as for the HETGS data, they are produced likely due to the use of the 
abbreviated data to model this feature. 

We note that some low ionization lines show residuals consistent with an outflow velocity 
of -1040 km s _1 as shown in the Figure [T51 This may indicate the presence of a HV-LIP (see 
§1]). However, attempting to model these residuals with an additional absorbing component 
does not improve the fit in a statistical way. 

4. Discussion 

The best fit model of the WA for the Seyfert 1 galaxy NGC 5548 requires two outflow 
velocity systems, each with two different ionization phases Figure [TBI shows that the 

profiles of the absorption lines (produced by ions such as Ne x, Fe xvill, O vm, and Mg xi) 
are fit better with two narrow line profiles than with a single broad component. Each of these 
two velocity components requires absorption by gas with very different ionization level. The 
LV component consists of gas with two different ionization states (the LV-LIP and the LV- 
HIP). The HV component requires also two absorbers with different ionization (the HV-SHIP 
and the HV-HIP). These four absorbers, divided into two velocity systems, can reproduce 
all the spectral features observed in the spectra (FeVll-xxiv, O vi-vm, Neix-x, Mgix-xil, 
Si XII- xiv, Cv-Vl, Six- XII, and Arxi-xn). We note that the HV system might also have 
an additional, low ionization component (HV-LIP). Figure [T5] shows a double line profile 
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for several low ionization level lines, suggesting a LIP with higher outflow velocity (-1041 
km s _1 ). While the inclusion of this component is not statistically required by the data, its 
physical properties appear to be well constrained (as can be observed in in Figure [THl where 
the 1, 2, and 3a confidence regions for the ionization parameter vs. the column density of 
this component are presented). Further evidence for the presence of this component arises 
from the fact that UV system 1 can be associated with the HV absorber (see §4.ip . and the 
HV-SHIP and HV-HIP have a very hgh ionization state to produce such absorption ((see 
Tables [7] and [S]). We conclude that the data clearly shows a shift to higher ionization material 
at higher velocity in NGC 5548 WA, with the LIP becoming weak and the SHIP dominating 
in the HV system. On the other hand, the LV system does not show the presence of the 
SHIP. 

The separation between components is a new and important feature in X-ray spectra of 
ionized absorbers, and could be achieved because of the large difference in velocity between 
the two systems (~ 500 — 600 km s" 1 ), as well as the different dominant ionization levels 
between them (see also Steenbrugge et al. 2005). 

4.1. UV Counterparts 

In the study of the narrow Ovi, Nv, and Civ absorption lines in the UV region, Cren- 
shaw et al. (2003) finds five separate kinematic components in HST -STIS and FUSE 
spectra. Our analysis reports two kinematic components (HV and LV systems). The HV 
system shows similar outflow velocity to the UV absorbing component 1 (—1041 km s _1 ). 
According to Crenshaw et al. (2003)b, the gas producing absorbing component 1 should pro- 
duce only a weak X-ray WA, which is consistent with the marginal detection of the HV-LIP 
in the Chandra data (§H]). 

There is also a correspondence between the X-ray LV system and UV absorbing compo- 
nents 2 and 3 (with outflow velocities —667 and —530 km s _1 , respectively). Crenshaw et al. 
(2003)b and Crenshaw et al. 2009 suggest that at least component 3 (and possibly compo- 
nent 2) should produce absorption in the X-ray range. This is consistent with the detection 
of the LV-LIP. There is no X-ray detection of UV absorbing systems 4 and 5 (with velocities 
—336 and —166 km s _1 respectively), although we cannot rule out a possible contribution 
from component 4 to the LV-LIP, given the uncertainty in the outflow velocity of the X-ray 
systems. According to the analysis of the UV data by Crenshaw et al., system 5 has the 
lowest column density (logN#=19.27). Krongold et al. (2003) discuss how the difference in 
sensitivity of the current X-ray and UV detectors may affect the detection or non-detection 
of a given UV component in the X-rays. Their Figure 11 shows the observational lower 
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limits for the equivalent widths of the CivAA1548.2, 1551, OviAA1038, 1032, O VIIA21.602, 
O VIIIA18.968 lines detectable in the UV and X-ray bands, as a function of the ionization 
parameter and H equivalent column density. Thus, given their ionization state and column 
density, systems 4 and 5 are consistent with not being detected in the X-rays, given the 
limited sensitivity of the X-ray observations. 

While a kinematic connection between the components is evident, an important question 
is how the predicted/measured column densities for the different analyses compare between 
them. The ionic column densities predicted by our model are now listed in Table [7] and [S] 
for the LV and HV systems, respectively. In these Tables, we also list the column densities 
predicted (or measured) by Crenshaw et al. (2009), based on UV data, and by Steenbrugge 
et al. (2005), based on the LETG 2002 data (their model B, see §4.31) . The column densities 
for ions CIV, NV y OVI between our predicted values and those measured by Crenshaw et al. 
(2009) in the UV differ by factors of 16, 4 and 1.2, respectively. We note that OVI (where 
no diference is found) is detected in the X-ray analysis as well. The differences between the 
UV and X-ray columns could be produced because of 3 different factors: 1) The blending 
of UV components 2 and 3 in the X-rays (and maybe an additional contribution from UV 
component in our LV-LIP model. 2) The different SEDs used in the 2 analysis: The SED 
used by Crenshaw et al. includes much more photons in the far (unobservable) UV, and 3) A 
detailed analysis on the dependence of the covering factor on velocity could not be performed 
in the analysis of the UV line profile. This could underestimate the columns as shown by 
Arav et al. 2003. Taking into account these effects, we consider that the predictions of 
our model for the columns measured in the UV are very reasonable, and conclude that the 
LV-LIP component can produce the absorption features produced by UV component 3 (and 
maybe 2 and 4). However, we note that the values predicted by Crenshaw et al. for the 
high ionization charge states differ from those measured in this analysis. The reason for this 
discrepancy, is that these ions are not observable in the UV, and only three ions (CIV, NV 
y OVI) were used in the UV model. 

With these considerations, and given the strong coincidence in kinematics between the 
UV and X-ray outflows, we conclude that a common origin is not only possible, but quite 
likely. Such connection was originally suggested by Mathur et al. (1994) for NGC 5548, and 
has also been suggested for several other objects (NGC 985- Krongold et al. 2005; NGC 
3783-Gabel et al. 2003, Krongold et al. 2003, Netzer et al. 2003; NGC 7469- Blustin et al. 
2003, Kriss et al. 2003). Thus, these studies on the WA may represent a general feature of 
the structure for quasars, as proposed by Elvis (2000). 
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4.2. Two Multi-phase Absorbing Outflows 

The presence of two different absorbing components with different temperature but 
similar outflow velocity suggests that the absorber may be formed by a multi-phase medium 
(e.g. Elvis, 2000; Krongold et al. 2003). Given the two velocity systems in NGC 5548, we 
could be seeing the presence of two different outflows (with different outflow velocity), each 
being a multi-phase medium. 

This is further supported by the pressure equilibrium of the two absorption components 
in each velocity system. In Figure [T7] we show the thermal stability curve of the gas (also 
known as the "S-curve", Krolik et al. 1981) for the SED used in our analysis (§3.2.1; Fig. 
[3]). The S-curve marks the points of thermal equilibrium in the T vs. U/T plane, where T is 
the photoionization equilibrium temperature of the gas, and U /T is inversely proportional 
to the gas pressure (U/T oc P ra d/ P gas )■ 

In Figure ITTa. the position on the S-curve of the two LV absorbing components is shown. 
The LV-HIP and the LV-LIP have quite different temperatures but, within the errors, are 
consistent with a single value of the U/T, and so can be in pressure equilibrium. Figure [TTb 
shows the HV-SHIP and HV-HIP on the S-curve. Again, the pressure equilibrium between 
the two components is allowed. Note that, given the possible range of ionization parameters 
of the tentative HV-LIP (Figure [TH]), this component, if indeed present (as suggested by our 
results and by those by Crenshaw et al. 2009), could also be in pressure equilibrium with 
the higher ionization ones that form this system. Thus, both the HV and the LV systems 
are formed by a two-phase medium (or maybe three-phase for the HV system), with one 
component pressure confining the others (see Krongold et al. 2005a for a more detailed 
explanation) . 

NGC 5548 thus joins the group of AGNs for which pressure equilibrium applies: NGC 
3783 (Krongold et al. 2003, Netzer et al. 2003), NGC 985 (Krongold et al. 2005, 2008), 
NGC 4051 (Krongold et al. 2007), Mrk 279 (Fields et al. 2007), UGC 11763 (Cardaci et al. 
2009). No other well studied WA in a Seyfert galaxy is inconsistent with the description of 
a multi-phase medium. We note however, that we will further test this idea for the WA over 
a large sample of objects (Andrade- Velazquez et al. in preparation). NGC 5548 is the first 
AGN in which a two-phase WA medium is likely for two velocity systems. The presence of a 
multi-phase medium obviously requires that the different ionization components that form 
it lie at the same distance from the ionizing continuum. We do not have the means to test 
this requirement (though see Krongold et al. 2010). However, for NGC 4051, Krongold et 
al. (2007) measure the distance of the two absorbing components and found that they were 
indeed consistent with being co-located, giving further plausibility to multi-phase absorbing 
winds. 
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The SHIP of the HV system and the HIP of both the LV and HV systems lie on unstable 
regions of the S-curve, implying a possible inconsistency with the multi-phase scenario. 
However, the shape of this curve depends on the exact shape of the SED of the source, 
which has major uncertainties, particularly in the extreme UV region (Haro-Corzo et al. 
2007). The chemical composition of the absorbing gas has an important effect on this curve 
also (see Komossa & Mathur, 2001, Fields et al. 2007). In addition, recent studies of the 
S-curve with new and more reliable dielectronic recombination rate coefficients (Badnell et 
al. 2006) show that there is a larger probability of having a thermally stable WA at 10 5 K, 
(Crenshaw et al. (2009) obtain a T ~ lxlO 5 for the UV component 3), and then a greater 
plausibility for its multiphase nature (Chakravorty et al. 2008). 



4.3. Comparison with Previous Works 

Steenbrugge et al. (2005) analyzed three Chandra grating observations of NGC 5548 
(those carried out on 2002, Table 1) with a total exposure time of 494 ks. They presented a 
model showing that a continuous distribution of column densities and ionization parameters 
could fit the data. They further concluded that the data was not consistent with discrete 
phases of material in pressure equilibrium. 

This conclusion is in contradiction to the results presented here, making a detailed 
comparison of both analysis necessary. In their model C, Steenbrugge et al. (2005) explored 
a model similar to the one presented here. They used a self-consistent model that takes 
advantage of photoionization balance calculations to make a global fit to the data. They 
concluded that following this approach, and fitting Fe separately, a good fit could be obtained, 
but constraining only two absorbing components. They further tested the presence of discrete 
components by adjusting equivalent H column densities and ionization parameters to the 
absorbing components of their model B, where ionic column densities were measured on an 
ion by ion fit. In this case they find that at least five different components with different 
ionization state were required. This approach led them to conclude that discrete components 
were not a good representation of the data. 

The most fundamental difference between the analysis carried out by Steenbrugge et 
al. (2005) and the present analysis is that their models B and C used only a single velocity 
outflow with a large turbulent broadening to fit the data. However, as they report in their 
paper, and as discussed here, the data are better fitted with two different velocity compo- 
nents, the HV system and the LV system. Furthermore, as shown here (and as also noted by 
Steenbrugge et al.) the degree of ionization is different for the different velocity components, 
as the HV system becomes dominant for high ionization states. This led them to conclude 
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that more and more ionization components were required to fit the data and thus a struc- 
ture of discrete clumps of gas was unlikely. However, here we have shown that, when each 
velocity component is modelled independently, both can be well described by two ionization 
components, making the presence of discrete components plausible and likely. Furthermore, 
the addition of extra components does not improve the fit. So, given the difference between 
number of free parameters from the different models, we can claim by Occam's razor that 
our global approach is preferable. Table [7] and [8] show the column densities for our model 
and the Steenbrugge et al. model B. The differences found between our analysis and that by 
Steenbrugge et al. (that analysis shows sistematically larger columns than ours) are likely 
produced because only one velocity component was used by Steenbrugge et al. to model the 
line profile produced by the two velocity systems. A detailed fitting including the line profile 
of each velocity system produces more reliable column densities. 

Steenbrugge et al. (2005) further tested whether their five discrete components could 
be in pressure balance, and they concluded that they could not. In fact, their S-curve (their 
Fig. 10) has no region where the slope is negative, making it impossible for the gas to be in 
pressure balance. We stress again that the shape of the S-curve depends strongly on the 
metallicity of the absorbing gas, and the shape of the SED used to produce it (in this work 
we use the best multiwavelength data available extracted from NED). Given this, and the 
new results by Chakravorty et al. (2008), we conclude that the possibility that the different 
phase are in pressure equilibrium cannot be discarded. Here we have shown that when the 
velocity systems are modelled independently, each one is consistent with pressure balance. 
We further note that Chelouche (2008) also find that the WA in NGC 5548 can be described 
by a multiphase medium. The author models the ionized flow assuming a thermal wind. 
While his model is more physically driven than ours, it does not take into account the two 
different velocity systems revealed by the data. However, given the similar results, the idea 
of a multi-phase medium absorber is reinforced. 



4.4. Possible Geometry of the Ionized Outflow 

The presence of two different velocity systems, each formed by gas with very different 
ionization state, can give a clue about the structure of the flow. In a continuous radial 
range of ionization structures, the wind must extend in our line of sight for parsecs (or 
more), forming a large scale outflow along our line of sight. We find difficult to reconcile 
such a distribution of material with the two velocity systems found in NGC 5548, even if 
inhomogeneities in the flow (like those suggested by Ramirez et al. 2005) are invoked. The 
presence of the two velocity components, then, suggests that we are not seeing the wind in 
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a radial direction. 

On the other hand, the presence of the different kinematic components can be naturally 
explained if we are looking at an outflow in a transverse direction, which is consistent with 
the detection of transverse motion in the UV absorbers (Marthur et al. 1995, Arav et al. 
2002, Crenshaw et al. 2003). In this case, an important component of the velocity (and 
acceleration) of the wind could take place in the direction perpendicular to our line of sight. 
Such a configuration for the wind is further supported by the finding that the different phases 
observed at each velocity form a single multi-phase outflow. 

If the wind arises from the inner accretion disk (Konigl & Kartje, 1994, Elvis et al. 
2000, and Krongold et al 2007), and flows in a transverse direction with respect to our 
line of sight, the two velocity components could be simply explained as forming in different 
regions of the disk (Elvis, 2000). The two velocity components in NGC 5548 are consistent 
with originate in the accretion disk, as discussed by Krongold et al. (2010). Furthermore, 
Crenshaw et al. (2009) suggest that UV component 3 (clearly associated to the LV-LIP), 
should be accelarated by magnetocentrifugal forces, and thus should form at an accretion- 
disk scale (Bottorff et al. 2000). We note that an accretion- disk wind does not imply a small 
scale outflow. The wind could still extend for large scales, but our line of sight would cross 
it only at a given distance from the central ionizing source. 

5. Conclusions 

We find that the WA of NGC 5548 can be modelled in a simple way with two different 
velocity systems (V out = -1110 km s _1 and V ou t — —490 ±49 km s -1 ) each with two phases of 
gas in pressure equilibrium. This results strongly suggests that the structure of the absorber 
in NGC 5548 could consist of two different outflows (with different outflow velocity), each 
formed by a multi-phase medium in pressure balance. The ionized absorber in NGC 5548 
now fits a pattern established for the other well studied systems. There are no remaining 
cases in which a multi-phase medium is inconsistent with the observations. 

The presence of two velocity systems, each formed by a multi-phase medium suggests 
that we are looking at the flows in a transverse direction with respect to our line of sight (as 
also suggested by the UV data). 

The velocity of the two X-ray absorbing systems is in agreement with components 1,2, 
and 3 found in the UV band (Crenshaw et al. 2003b). The X-ray data are further consistent 
with the non detection of UV components 4 and 5, given the limited sensitivity and resolution 
of the current X-ray (compared to those in the UV). There is also a reasonable agreement 
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in the column densities of CIV, NV, and OVI predicted by our X-ray analysis, and those 
measured in the UV. Thus, our models elegantly solve the long-standing problem of apparent 
discrepancies between the UV and X-ray absorbers, and give further support to the idea that 
UV and X-ray absorbers are part of the same outflow. 

This work was supported by the UNAM PAPIIT grant IN 11 8905 and the CONACyT 
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Table 1. Observations Log of NGC 5548. 



HETG/ACIS-S 



Observation 


UT Start Date 


Exposure 




S/N(22A) 


(2)837 


2000-02-05 15:37 


82.3 


1.14 


6.63 


(3)3046 


2002-01-16 06:12 


153.9 


1.96 


4.69 


TOTAL 




236.2 




8.12 


LETG/HRC-S 


(1)330 


1999-12-11 22:51 


85.98 


2.99 


34.27 


(4)3045 


2002-01-18 15:57 


169.68 


1.38 


35.71 


(5)3383 


2002-01-21 07:33 


171.02 


2.19 


33.32 


(6)5598 


2005-04-15 05:18 


116.43 


0.56 


14.45 


(7)6268 


2005-04-18 00:31 


25.19 


0.46 


9.38 


TOTAL 




568.3 




62.09 



a In ks 

b In 10 _11 erg s _1 cm -2 



Table 2. Models used to fit the NGC 5548 spectra. 



M 


Power Law 




Black Body 




EL d 




WA 


* 




Goodness 




r 


A , a 
^-pwlw 




kT (keV) 


A BB b 






SHIP^y 




HIPiv 




X\ dof, ^ 


A 


1.52±0.006 


4.87±0.02 x 10" 


-3 


0.012± c 


1.00xl0- 9 ± c 




4 










3813, 4112, 0.93 


B 


1.560±0.007 


5.20±0.03 x 10- 


-3 


0.085±0.005 


6.98±1.36 x 10- 


-5 


4 




V 






3053, 4108, 0.74 


C 


1.568±0.007 


5.00±0.03 x 10" 


-3 


0.083±0.005 


7.28±1.54 x 10" 


-5 


4 


V 


V 






2987, 4104, 0.73 


D 


1.595±0.008 


5.00±0.05 x 10- 


-3 


0.099±0.006 


5.65±1.07 x 10- 


-5 


1 


V 


V 


V 




2946, 4100, 0.72 


E 


1.597±0.009 


5.40±0.05 x 10- 


-3 


0.110±0.007 


5.17±0.85 x 10- 


-5 


4 


v' 


V 


v' 


v' 


2904, 4096, 0.71 


F e 


1.598±0.010 


5.46±0.05 x 10- 


-3 


0.116±0.009 


4.51±0.75 x 10- 


■5 


8 


v' 




v' 


V 


2885, 4102, 0.70 


G 


1.778±0.02 


7.84±0.05 x 10" 


-3 


0.083±0.002 


5.93±0.02 x 10" 


■5 


7 


V 


V 


V 


V 


2013, 1618, 1.24 



a In units of 10 -3 ph s _1 cm -2 keV -1 at 1 keV 

b In units of 10 — 5 Lsg/Df , where L39 is the source luminosity in units of 10 39 erg s _1 and D±o is the distance to the source in units of 10 kpc 
c The uncertainties cannot be found due to the presence of the absorber. 
d EL means Emission Lines. 

° Model E plus the other lines remains (see Table [3j|. 
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Table 3. Emission lines in the NGC 5548 spectrum 1 



ION Rest Frame A(A) b Measure A(A) FWHM (kms^ 1 ) EW(mA) 



HETG SPECTRUM 



Fe K Q 


1.937 


1 Q3fiS +a0086 
±.youo_ 0070 


2644.4±2469 


27.9 ttli 


Si xiii 


6.648 


6.739iS:88g 


754±938 


q 7 +1.4 
' -1.4 


Ne ix 


13.699 


13.678 ±§:88| 


485±593 


22.6± 317 


O vn 


21.607 


Z.L.OJI _ .oi8 


936±636 


oo.y_ 19 7 


O vn 


21.807 


21 785 +a005 

AL.IOQ -0.007 


928±299 


73.2±H:| 


O vn 


22.101 


22 075 +a007 

ZZ.UIO -Q.008 


915±311 


244.2 + |ti 


Nvi 


23.024 


23.035 i8$S 


327±373 


26.21^:8 


Nvi 


23.277 


23.306i8:8}8 


673±398 


80.91^ 



LETG SPECTRUM 



Ne IX 


13.699 


13.677 18'° 16 


60.7±822 


5.9± 2 .2 


O vn 


21.607 


21 623 +a022 

Zl.UZO _o.019 


944±513 


34.9 + 7 '| 


O vn 


21.807 


21 760 +a025 

Zl. (UU -Q.013 


935±554 


42.8l 23 6 1 


O vn 


22.101 


22 075 +ao ° 4 

ZZ.LWd -0.004 


923±159 


ii9.o±ig:i 


Nvi 


29.534 


29 502+ 021 


462±884 


23.9™ 


Nvi 


29.081 


29.029 c 


10 c 


H.Slo 7 - 1 


Cv 


41.472 


41.429 c 


20±452 





a All line positions are reported with respect to the rest frame of the host 
galaxy. Given the uncertainties in the measurement, as well as the filling of 
the absorption lines by the emission features, all lines are consistent with 
the rest frame of the host galaxy. 

b The line rest wavelengths were taken from the NIST Atomic Spectra 
Database (URL: http:/ /physics. nist.gov/PhysRefData/ASD/index.html) 

c The uncertainties cannot be found. 
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Table 4. The best continuum parameters. 



Power Law 



Photon Index (r) Normalization 11 N^Gai (cm 2 ) 
HETGS 1.59±0.010 54±0.5 1.83 x 10 20 

LETGS 1.77±0.02 78±0.5 



Thermal component 

kT (keV) Normalization 13 
HETG 0.11±0.009 0.45±0.08 
LETG 0.083±0.0016 0.59±0.002 



a In 10 4 photons keV 1 cm 2 s 1 at 1 keV. 

b In 10~ 4 L 39 /Df , where L 39 is the source luminosity 
in units of 10 39 erg s _1 and D\q is the distance to the 
source in units of 10 kpc. 



Table 5. Models B-D Parameters for the Ionized Absorber 



COMPONENT 


logU 


logNjj(cm 2 ) 


Vel out (km s x ) 


Vel tur6 (km s l ) 


Bl 


0.75±0.02 


21.58±0.03 


-740±150 


536±35 


CI 


1.06±0.10 


21.59±0.06 


-1120±150 


117±26 


C2 


0.63±0.02 


21.28± 0.04 


-450±150 


241±37 


Dl 


1.14±0.02 


21.71±0.08 


-1120±150 


113±42 


D2 


0.67±0.02 


21.30±0.04 


-450±150 


275±117 


D3 


-0.49±0.08 


20.74±0.10 


-590±150 


105± 47 



Table 6. Final Best Fit Ionized Absorber Parameters (Model E and F). 



Parameter V/ l j=— 1110 km s 1 V; =— 490 km s 1 





SHIP 


HIP 


HIP 


LIP 


Log U a 


1.23±0.06 


0.67±0.03 


0.67±0.03 


-0.49±0.09 


Log N H (cm" 2 )" 


21.73±0.12 


21.03±0.07 


21.26±.04 


20.75±0.10 


V T urb (km s" 1 ) 


175 ±67.0 


100 ±23.0 


177 ±39.0 


105 ±29.0 


Vout (km s- 1 )" 


-1040±150 


-1180±150 


-400±150 


-590±150 


T (K) 6 


27.0 ±5.6 


5.8 ± 1.0 


5.8 ±0.8 


0.35 ±0.04 


[Log T (K)] 


6.47±0.12 


5.89±0.06 


5.72±0.12 


4.54±0.02 


Log T/U (ocP c ) 


5.21±0.17 


5.15±0.15 


5.04±0.14 


5.06±0.1 



a Free parameters of the model. 

b In units of 10 5 Kelvin. Derived from the column density and ioniza- 
tion parameter, assuming photoionization equilibrium. 

°The pressure Pocn e T. Assuming that both phases lie at the same 
distance from the central source n e ocl/U, and PocT/U. 
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Table 7. Ionic Column Densities for the LV system 



ION 


logNzjp 


logN^/p 


\ogNsuM 


logN iC 


logN i5 






IDT'? 








r v 

v_> V 


lUi 1 o 




1 7^ 

-LU . I o 




1 7 90 


c vt 

V_. V 1 


1 7 0^ 

X 1 .WO 


1 D3 


1 7 (17 

X 1 .U I 


is m 

i o. iu 


1 7 fin 

1 l . UU 


Nv 

IN V 


1 ^ ^ 


1 1 07 

ll.U i 


1 ^ ^ 

10.00 


14 cn 

11. c/O 




N VT 
IN VI 


1 R 4R 
XU.^U 


1 4 Oft 


1 4fi 
1U.1U 


1 7 90 


1 7 90 
1 ( .zu 


N VTT 
IN V 11 


1 R 9R 


1 ^ Rfi 
1 u. ou 


1fi 41 


1 R 


1 7 ^10 

1 I .OU 


O v 

V J V 


1 R ^R 


1111 

1 1 . -1 1 




1 ^ 70 


1 R0 
1 u. ou 


OVT 
V 7 V 1 


1 R KQ 

XU. Ot/ 


1^1^ 
10. 10 


1 RQ 


1 R0 


1 fiO 
1 u. uu 


O VTT 


1 7 40 

X f i^lU 


1 ^ R& 


1 7 41 

1 / .ft 1 


1 R ^0 


1R 1R 


W V 111 


XU. i O 


1 7 OR 


1 7 94 


1 R 70 

lO. i u 


1 R ^ 


Np TV 

IN t: 1 A. 


1 ^ Q9 


1 R 1 


Ifi ^9 


1 R 00 
1 o. uu 


1 7 90 


Np V 




1 R 

X U iOO 


1 8R 


1 7 fiO 
1 i .uu 


1 7 Q0 


Mgxi 


15.04 


16.20 


16.23 


17.30 


17.00 


Mg xii 


13.31 


16.51 


16.51 


16.50 


17.30 


Si vm 


15.46 


12.06 


15.46 


16.70 


16.20 


Si ix 


15.83 


13.40 


15.82 


17.20 


16.00 


Si xi 


15.44 


14.34 


15.47 


16.80 


16.00 


Si xm 


13.81 


16.45 


16.45 


16.40 


17.00 


Si xiv 


11.74 


16.39 


16.39 


15.20 


17.00 


Sxi 


15.35 


14.99 


15.51 


16.90 


16.00 


Sxn 


14.84 


15.63 


15.70 


16.70 


16.20 


Fe xvn 


13.82 


15.73 


15.74 


15.10 


16.30 


Fe xviii 


13.54 


16.19 


16.19 







a Columns 2 and 3 correspond to the ionic column densities 
predicted by the LV system of the best model in this work. 
Column 4 is the added contribution of the two ionization 
components (LV-LIP and LV-HIP) of this system. The pre- 
dicted ionic column densities by UV component 3 (Crenshaw 
et al. 2009) are located in column 5. Column 6 shows the 
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column densities reported by model B of Steenbrugge et al. 
2005. All columns are in units of cm -2 . 
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Table 8. Ionic Column Densities for the HV system 



ION 


\ogN HIP 


logNsHIP 


logN sum 


logN iS 


Civ 


10.36 


9.04 


10.38 





Cv 


12.98 


12.74 


13.18 


17.20 


Cvi 


16.26 


15.63 


16.35 


17.60 


Nv 


11.30 


9.57 


11.31 




Nvi 


14.30 


12.96 


14.32 


17.20 


Nvn 


16.09 


15.50 


16.19 


17.50 


Ov 


10.86 


9.30 


10.87 


16.80 


Ovi 


13.37 


11.30 


13.37 


16.60 


vii 


15.90 


14.51 


15.92 


18.18 


vm 


17.29 


16.76 


17.40 


18.53 


Ne IX 


15.86 


14.68 


15.89 


17.20 


Nex 


16.65 


16.55 


16.90 


17.90 


Mgxi 


15.96 


15.10 


16.02 


17.00 


A/Trr VTT 
IVlg, .All 




1 AA 


1U.U i 


1 i .ou 


Si vm 


11.83 


7.98 


11.83 


16.20 


Si ix 


13.16 


7.98 


13.16 


16.00 


Si xi 


14.10 


11.85 


14.10 


16.00 


Si xin 


16.22 


15.75 


16.35 


17.00 


Si xiv 


16.16 


16.71 


16.82 


17.00 


Sxi 


14.17 


9.71 


14.17 


16.00 


Sxn 


14.75 


12.37 


14.75 


16.20 


Fe xvn 


15.49 


13.24 


15.49 


16.30 


Fe xvin 


15.96 


14.12 


15.97 





a Columns 2 and 3 correspond to the ionic column 
densities predicted by the HV system of the best 
model in this work. Column 4 is the added contribu- 
tion of the two ionization components (HV-HIP and 
HV-SHIP) of this system. Column 6 shows the col- 
umn densities reported by model B of Steenbrugge et 
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2005. All columns are in units of 
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Fig. 1. — NGC 5548 light curve for the observations analyzed in this work. Points 2 and 3 
correspond to HETGS observations and points 1, and 4-7 to LETG observations. 
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Fig. 2.— Model A plotted over Chandra-MEG data ( [1.6-24.5] A) of NGC 5548. The model 
(red thick line) consists of a powerlaw and a black-body component plus four strong emission 
lines (blue labels). Possible transitions due to the ionized absorber are indicated with red 
labels. 
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Fig. 3. — Spectral Energy Distribution used in this analysis (black solid line) to model the 
ionized absorber in the Chandra data of NGC 5548 (the SED is based on our observations 
and NED data, see §3] for details). Blue labels mark different observational bands (see § 
3). The slopes adopted in each band are marked with red labels (this slope relates to the 
photon index as T=s+2). The dashed green line represents the SED used by Steenbrugge et 
al. 2005. 
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Fig. 4.— Model B plotted over Chandra-MEG data ( [1.6-25] A). This model (red thick 
line) includes the same components of model A plus one absorbing component. Emission 
and absorption lines are marked as in Figure HJ 
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Fig. 5. — Model B and residuals in selected regions of the spectrum showing that the presence 
of additional absorbing material with higher and lower ionization is required by the data. 
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Fig. 6. — Detailed line profile for several absorbing features that are not well fitted by only 
one broad absorbing component. The dotted lines correspond to the outflow velocity of the 
UV components (Crenshaw et al. 2003). Significant residuals are found with outflow velocity 
consistent with UV component 1. 





Fig. 7. — As Figure HI but for model C, which includes two absorbing components. 
This model can fit transitions by both high and medium ionization level ions such as 
Fexxni-xxiv- 
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Fig. 8. — As Figure HJ but for model D (red thick line), which includes three absorbing 
components. The contribution of the lower ionization component is observed in transitions 
such as Mgix-xi and Ovii- 
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Fig. 9. — Model F (red thick line; corresponding to model E plus 4 additional emission 
lines, see § 3.4 ) plotted over the MEG spectra ([5-24.5] A) of NGC 5548. Two velocity 
components, each formed by two phases, are required to fit the data. The upper labels 
correspond to the absorption lines, while the lower labels mark the emission features. 
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Fig. 10.— Model F (red thick line) plotted over the HEG spectra ([1.6-8.84] A) of NGC 5548. 
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Fig. 11. — MEG spectrum and model between (5 and 24.5) A showing only the contribution 
of the low outflow system. The left panel shows the contribution by the Higher Ionization 
Phase (HIP) and the right panel shows the Low Ionization Phase (LIP). 
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Fig. 12. — MEG spectrum and model between (5-24.5) A showing only the contribution of 
the HV-HIP, and the right panel the HV-SHIP. 
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Fig. 13. — Model F (blue thick line) plotted against selected HEG and MEG spectral regions, 
showing the best fit lines profiles. The two outflow velocity systems are marked with letter 
E at the top of the panels, and correspond to -1180 and -400 km s _1 . The fit with a single 
broad profile (model B, red thick line) is also shown. The corresponding velocity of this 
model (-740 km s -1 ) is marked with the letter B and a dashed vertical red line. The UV 
velocity systems (Crenshaw et al. 2003) are marked with dotted lines. 
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Fig. 14. — The best fit for WA on the HETGS data extrapolated to the LETG data (model 
G). This model shows a good fit on the LEG spectral region ( [9-41.6] A). The O vn , N VI , 
Ne/x and Cy emission lines are marked with dotted blue lines. 
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Fig. 15. — As Figure El but over the LETGS spectrum. The line profiles of the low ionization 
phase suggest the presence of 5th absorbing component, i.e. a low ionization phase at high 
outflow velocity. 
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Fig. 16. — Confidence regions (at 1, 2, and 3a) for the ionization parameter vs. column 
density diagram constrain the possible values of these for High Velocity - Low Ionization 
Phase (logU=-0.47 and logN#=20.35). 
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Fig. 17. — Thermal equilibrium (S-curve) for the SED used in the present analysis. Left: 
The two absorbing components forming the Low Velocity system are plotted in the log(U/T) 
vs log(T) diagram, where log(U/T) oc log(P~ 1 ). Right: The diagram for two absorbing 
components forming the High Velocity system. In both systems the pressure balance between 
the components is evident. The values for HV-LIP might be consistent with pressure balance 
as well. 



